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The origin of the correlation between the UV and X-rays 
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ABSTRACT 

We propose that the UV/X-ray correlation observed in NGC 4151 by EXOSAT and IUE 
is due to reemission of the X-ray flux absorbed by cold clouds in the line of sight. This 
transmission model satisfies the energy balance, provides a good fit to the X-ray and 
UV data, and predicts no Compton reflection of X-rays (as observed) . The correlation 
can alternatively be explained by reprocessing of X-rays and 7-rays by a cold accretion 
disk. A specific model in which X-rays and 7-rays are emitted by a dissipative patchy 
corona and have the spectrum cut off at ~ 100 keV (as seen by OSSE) is shown to be 
energetically possible even if it does not give as satisfactory a fit as the transmission 
model. However, the model predicts a Compton- reflection spectral component, whose 
presence remains to be tested by future detectors sensitive above 10 keV. 
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1 INTRODUCTION AND SUMMARY 

Correlation between the UV and X-ray fluxes in radio-quiet 
AGNs ranges from good to non-existent (e.g., Ulrich 1994). 
So far, a good correlation was observed in the Seyfert galax- 
ies NGC 4151 (Perola et al. 1986, hereafter P86) and NGC 
5548 (Gavel et al. 1992). Here we consider the origin of the 
correlation in NGC 4151, a nearby Seyfert 1.5 galaxy at the 
distance of D ~ 20 Mpc. 

The correlation in NGC 4151 was observed during two 
IUE/ EXOSAT campaigns in 1983 November 7-19 and 1984 
December 16-1985 January 2 (P86). The UV flux at 8.5 
eV was proportional to the 2-10 keV X-ray flux corrected 
for absorption (Section^). The time delay between the flux 
variations was less than ~ 2 days. The delay, however, was 
more tightly constrained to 5s 0.5 day in 1993 December 
(Warwick et al. 1996), when the light curves were sampled 
at time intervals much shorter than in 1983-86. 

In the first part of this Letter, we present a new model 
explaining the correlation. We propose that the UV radia- 
tion is from reemission of X-rays absorbed by cold matter 
outside the X-ray source (Section [|). Thus, we identify the 
matter responsible for the UV emission with that responsible 
for the observed absorption of X-rays. The column density, 
Nn ~ 10 23 cm -2 , is determined from modeling the X-ray 
spectra (P86). We show that this model (hereafter referred 
to as the transmission model) gives an excellent fit to the 



UV and X-ray data as well as it satisfies the energy balance. 
The model implies only residual Compton reflection (Light- 
man & White 1988) because the absorber is Thomson-thin, 
r ~ 0.1, which is in agreement with hard X-ray observations 
(Maisack & Yaqoob 1991; Yaqoob et al. 1993; Zdziarski, 
Johnson & Magdziarz 1996, hereafter Z96). The absorbed 
flux is mostly from the observed 2-10 keV range and thus the 
model does not depend on the shape of the 7-ray spectrum, 
which was not observed during the X-ray/UV observations. 

From considering the absorber column density, ioniza- 
tion balance, and free-free absorption we determine that the 
absorber is not uniform but rather made of many dense dis- 
crete clouds. The clouds are optically thick to UV radiation 
and their temperature is kT ~ 3 eV. The obtained cloud pa- 
rameters are similar to those postulated to exist in the inner 
regions of AGNs by Guilbert & Rees (1988) and Ferland & 
Rees (1988). 

After presenting the transmission model we reexamine 
a previously proposed model of reprocessing by an optically- 
thick accretion disk (hereafter referred to as the reflection 
model; Section Q). In that model, the correlation is due to the 
disk irradiated from above by X-rays and 7-rays (hereafter 
X-7) and reemitting the absorbed flux in the UV (Ulrich et 
al. 1991; Ulrich 1994; Perola & Piro 1994, hereafter PP94). 

A number of problems appear in the reflection model. 
One is the observed absence of a strong spectral component 
at <; 10 keV expected from Compton reflection of the irra- 
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diating X-rays from the optically thick disk (Yaqoob et al. 
1993; Z96). [This contrasts the other Seyfert with a UV/X- 
ray correlation, NGC 5548, which X-ray spectrum does have 
a prominent Compton-reflection component (Nandra et al. 
1991).] Furthermore, the most detailed version of the model, 
proposed by PP94, requires that in order to satisfy the en- 
ergy balance the extrapolated X-ray emission has to extend 
to ~ 4 MeV, which has not been confirmed by recent OSSE 
measurements (Z96). In that case, the observed large intrin- 
sic dispersion of the 2-10 keV spectral index (P86) implies 
that the extrapolated power-law flux near ~ 4 MeV (which 
dominates the absorbed flux for the hard spectrum of NGC 
4151) becomes a strongly dispersed function of the 2-10 keV 
flux, which in turn implies a poor X-ray/UV correlation (see 
Section^), contrary to the data. 

We show here that most of the above difficulties can be 
dealt with by choosing a source geometry different from the 
central optically-thick spherical emission region proposed by 
PP94. Namely, we propose that the X-7 source forms a dis- 
sipative corona above the surface of the disk, as in the model 
of Seyfert l's of Haardt & Maraschi (1993) and Haardt, 
Maraschi & Ghisellini (1994). However, the corona is more 
patchy than in Seyfert l's to account for X-ray spectra of 
NGC 4151 being much harder than those of average Seyfert 
l's (Z86). We find that the corona geometry results in a 
strong reduction of the required X-7 luminosity with respect 
to the model of PP94. We find that the UV/X-7 energy bal- 
ance is satisfied even if the incident X-7 spectrum is cut off 
at ~ 100 keV, as observed by OSSE. This low cutoff energy 
also reduces the dispersion in the correlation of the total ab- 
sorbed power with the 2-10 keV flux. However, the model 
does predict a Compton-reflection spectral component only 
marginally consistent with hard X-ray observations. 

Note that both models are similar in that both pos- 
tulate absorption of X-rays by some cold medium, which 
subsequently reprocesses and reemits the absorbed flux as 
UV photons. The principal difference is that the medium is 
Thomson-thin in the transmission model whereas it is very 
Thomson-thick in the reflection model. 
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Figure 1. The correlation between the absorption-corrected EFe 
at 5 keV and at 8.5 eV observed during the 1983-1986 campaigns 
by EXOSAT and WE. The vertical error bars on Figs. 1, 2, 3a 
are for A\ 2 = 2.7, i.e., 90 per cent confidence for one parameter 
of interest. 



Fig. shows the resulting correlation between the 
absorption-corrected EFe at 5 keV (of the power-law only) 
and at 8.5 eV [from IUE, dereddened with E(B-V) = 0.05, 
PP94; Kriss et al. (1995) obtain a similar value of 0.04]. We 
see the correlation is close to linear at high significance. The 
values of EFe in X-rays are similar to those in the UV, 
which is indeed highly suggestive for the origin of the UV 
from some form of reprocessing of X-rays. 



3 THE TRANSMISSION MODEL 



2 THE DATA 

We have refitted the EXOSAT ME argon data (from 1.2 
to 10-15 keV) using a model consisting of a power law, 
Fe Ka line, absorption by a dual neutral absorber (i.e., a 
product of partial and complete covering, P86), and a con- 
stant soft X-ray component modeled as in Z96. We assume 
an Fe overabundance of 2.6, which is the average value ob- 
tained by P86. We find the X-ray energy spectral indices, 
ax {Fe oc E~ ax , where Fe is the energy flux), are weakly 
constrained by the argon data alone. Those indices were de- 
termined before in the 3.5-25 keV range by P86, who used 
the ME xenon data extending up to 25 keV in addition to 
the argon data. The addition of the xenon data resulted in 
the allowed ranges of ax much narrower than ours. Since the 
xenon ME data are no longer available we used the ranges 
of ax as given in P86 (in the range of ~ 0.2-0.7), which 
we found consistent with the argon data. (We have checked 
that relaxing this assumption increases the error bars shown 
on the figures here but it does not affect our qualitative con- 
clusions.) 



The fits in Section g give us the fluxes absorbed by the 
cold matter in the line of sight, F" a b s . We assume that the 
absorbing clouds are dense enough for almost complete ther- 
malization of that flux (see below). We determine the tem- 
perature of the absorber by requiring that the spectral in- 
dex in the UV between 8.5 and 7.2 eV equals the average 
observed value, a = -0.15 (Ulrich et al. 1991; PP94). This 
gives kT ~ 3 eV. Fig. |i| compares the predictions of that 
model with the data, assuming 



Fvv = fF ahE + F , 



(1) 



where Ftjv is the integrated UV flux (related to the plotted 
flux at 8.5 eV by the blackbody formula), Fo = 0.05 keV 
cm -2 s _1 is a residual flux, and / = 0.60 takes into account 
a covering of the X-7 source by the clouds of < 47r as well as 
an efficiency (< 1) of the conversion of the absorbed flux into 
the blackbody continuum. We see that there is a very good 
agreement between the data and the model in spite of our 
rather simplified assumptions (constant T, /, Fo, pure black- 
body UV spectrum). Since / < 1, the model does satisfy the 
energy balance, i.e., there is more than enough power in the 
absorbed X-rays to account for the observed UV emission. 
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Figure 2. Comparison of the prediction of the transmission 
model with the IUE/EXOSAT data. The crosses give the ab- 
sorbed flux, -F a t, s , obtained from the EXOSAT data, and EFe (8.5 
eV) is from IUE. The solid line gives EFe(8.5 eV) as a function 
of -F a bs from the transmission model (see text). 



We consider now constraints on the absorbing medium. 
For simplicity, we assume that all the absorption occurs in 
a single cloud. From the fits, we obtain the typical column 
density of a cloud in the partial absorber of Nn ~ 10 23 
cm -2 with the typical covering factor of ~ 0.5. The Nn of 
the complete absorber is a few times less. The dominant 
thermalization process at 8.5 eV is free-free, which implies 
that the size of a single cloud, r c , is Si 10 7 cm (equivalently 
n r& 10 16 cm -3 ), in order for the free-free optical depth 
at energies Si 10 eV to be > 1 at T and Nn as above. At 
those conditions, collisions alone are sufficient to almost fully 
ionize hydrogen, and thus there are enough free electrons for 
free-free absorption. 

The average size of the entire absorber, r, can be de- 
termined from the blackbody law, Luv (~ 10 43 erg s _1 ) 
= 47rtTV, to be ~ 10 14 cm, which satisfies the limit 
r SS 10 15 cm (for geometry close to spherical) from the rel- 
ative UV/X-ray time delay (Section |lj). The ionization pa- 
rameter, L/(nr 2 ), is then Ss 0.1, which implies metals are 
almost neutral, and thus soft X-rays are strongly absorbed 
in bound-free processes. The parameters of the absorber are 
similar to those studied by Ferland & Rees (1988; see their 
Fig. 8) and Guilbert & Rees (1988), for which they find the 
UV spectrum at Si 10 eV is indeed thermalized. Equiparti- 
tion magnetic fields can confine the clouds with the above 
parameters (Celotti, Fabian & Rees 1992). The X-7 source 
has to cover ~ 0.1 of the sky as seen from a typical cloud in 
order to satisfy the X-ray constraints on soft-photon starva- 



tion (Z96). Summarizing, the transmission model both ex- 
plains the data and appears to be self-consistent. 

It appears that accretion of the clouds in the observed 
absorber alone is insufficient to power the AGN. The total 
mass in the clouds is independent of their number, and is 
~ /(47r/3)A r H r 2 m p ~ 10 28 g. On the other hand, the mass 
required to account for L ~ (3-10) x 10 43 erg s _1 (assuming 
a cutoff at ~ 100 keV as observed by OSSE) at an efficiency 
of 0.1 is ~ (3-10) x 10 27 //3 r g, where (3 r is the radial veloc- 
ity. Thus, the two masses marginally agree only for highly 
relativistic radial velocities. 



4 THE REFLECTION MODEL 

Here we reexamine the reflection model, in which the correla- 
tion arises from X-7 irradiating an optically-thick accretion 
disk (e.g., PP94). The incident X-7 flux is mostly absorbed 
since the albedo is A ~ 0.3 (Z96). The absorbed power is 
reprocessed and reemitted locally by the disk mostly as UV 
blackbody radiation, similarly as in the transmission model. 

PP94 have developed a detailed model explaining the 
correlation, in which the X-7 source forms an optically thick 
central sphere irradiating a surrounding disk. They find that 
the amplitude of the UV emission at 8.5 eV and the 7.2-8.5 
eV spectral index (Ulrich et al. 1991) can be both explained 
only when the hard power-law emission extends to rather 
hard 7-rays; e.g., a cutoff above ~ 4 MeV is required for a 
typical q ~ 0.5. (We have reproduced those results using 
the assumptions of PP94.) 

As discussed in Section [l], there are major difficulties 
with such hard 7-ray emission. First, it has not been seen 
by the modern 7-ray detectors, OSSE and SIGMA (e.g., Z96; 
Finoguenov et al. 1995). Although 7-rays were not measured 
during the EXOSAT/ IUE campaign, a UV/X-ray correla- 
tion was also seen in 1993 December (see Section |E|; War- 
wick et al. 1996), when the OSSE data showed kT ~ 50±^ 
keV (Z96). Second, q x varies from 0.2 to 0.7 in the EX- 
OSAT data (P86). As a result of the dispersion of ax, the 
total fluxes, Fx-,, in the power-law spectra extrapolated to 4 
MeV correlate very poorly with the 2-10 keV fluxes. We find 
that indeed the corresponding reprocessed fluxes are away 
by factors up to 3 from their best-fit linear correlation with 
the UV, which invalidates the model (unless the dispersion 
is compensated for by suitable shapes of the 7-ray spectra, 
which seems unlikely). 

However, it is possible that other source geometries 
could be more compatible with the data. We study here 
a model in which the X-7 source forms a patchy corona 
above the surface of an accretion disk (Haardt, Maraschi 
& Ghisellini 1994). We assume that all dissipation takes 
place in the corona (Haardt & Maraschi 1993; Svensson & 
Zdziarski 1994). Thus, the local power, Q(r), released in the 
hot corona depends on radius, r, as the disk dissipation rate 
(Shakura & Sunyaev 1973), 
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/3rs\ 1/2 



(2) 



where r$ = 2GM/c? is the Schwarzschild radius, M is the 
black hole mass, and M is the mass accretion rate. The 
power released in the corona is emitted partly outward and 
partly towards the disk. We assume here that the luminosity 



4 A. A. Zdziarski and P. Magdziarz 



intercepted by the disk equals R times the luminosity emit- 
ted outward (and observed). R can be < 1 due to, e.g., a 
finite size of the cold disk or an outflow of the hot plasma. On 
the other hand, the luminosity emitted outward is related 
to the total observed X-7 flux, Fk 7 , by Lx-y = 47rD 2 Fx 7 . 
Then the luminosity absorbed and reemitted thermally by 
the disk is, 



(1 - A)RLx~, = 4tt(1 - A)RD 2 Fx~, 



(3) 



where A (~ 0.3) is the albedo. Energy balance then implies 
that the luminosity Luv equals the disk emission at the rate, 



a[T(r)] 4 = Q(r)(l-A)R/(l + R), 



(4) 



integrated over the disk surface. This allows us to solve 
for the constant, MM, in equation (^), as well as it shows 
that the UV spectrum is a standard blackbody disk spec- 
trum (Shakura & Sunyaev 1973) (rather than the single- 
temperature blackbody spectrum in the transmission model, 
or the spectrum from irradiation by an optically thick sphere 
in the model of PP94). 

The total observed UV flux is given by, 



Fuv = 2/i(l - A)RFxy, 



(•») 



where the factor 2/x accounts for disk emission with a con- 
stant specific intensity. Using the blackbody disk spectrum, 
we can then relate the observed 8.5 eV flux to the total UV 
flux, Fuv, as a function of r|/x, where fj, = cosi and i is the 
disk viewing angle (see PP94). Fig. shows this relation for 
several values of rs assuming /1 = 0.42 [as obtained by HST, 
Evans et al. (1993); rs for other values of scales as fi~ 1 ' 2 ]. 
Fig. shows the corresponding 7.2-8.5 eV spectral indices, 
compared with the observed values (Ulrich et al. 1991). 

The curves in Fig. can be compared with the to- 
tal fluxes, Fuv, absorbed by the disk and reemitted in the 
UV as obtained extrapolating the EXOSAT power laws to 
7-rays in the way observed by OSSE, i.e., for thermal Comp- 
tonization at kT = 60 keV (Z96). Z96 obtained the best fit 
values of R ~ 0-0.3 (at fi = 0.42) based on measuring a 
Compton reflection component in hard X-ray spectra, with 
R = 0.5 marginally consistent with the data (see also Mai- 
sack & Yaqoob 1991), and we thus use R — 0.5. Note that 
i? < 1 would not be consistent with the assumption that 
reprocessing is in the disk. Crosses in Fig. give Fuv for 
Rfi — 0.21 (corresponding to R = 0.5 and fx = 0.42). 

We can now compare predictions of the disk model for 
various rs with the UV spectral indices and Fuv inferred 
from the EXOSAT data. We find first that r s < 2.4 x 10 12 
cm fit the UV spectral indices, with the upper limit corre- 
sponding to Ax 2 = 2.7. The value of r s = 1.3 x 10 12 cm 
provides the best fit to the absorbed fluxes inferred from 
the X-ray data for R/j, — 0.21. The upper limit on rs corre- 
sponds to Rfi — 0.15. We note that the relative contribution 
of Compton reflection to the hard X-ray spectra decreases 
slower than oc fi (Magdziarz & Zdziarski 1995), and thus 
the problem of the weakness of Compton reflection together 
with the presence of reprocessing in NGC 4151 cannot be 
resolved by postulating < 1 (Z96). 

Note that the fit in Fig. |3p is significantly worse than 
that for the transmission model (Fig. ^). This is caused by 
the main contribution to Fuv being from ~ 100 keV, rela- 
tively far away from the 2-10 keV range, where the linear 
correlation is observed. Still, the correlation is much closer to 
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Figure 3. Comparison of the predictions of our reflection model 
with the data, (a) Crosses give the total U V flux, -Fuv , from ree- 
mission of the absorbed X-7 flux as obtained extrapolating the 
EXOSAT X-ray power laws to 7-rays assuming thermal Comp- 
tonization at kT = 60 keV, see text. The curves relate the total 
Fuv to FFg(8.5eV), as predicted by the disk spectrum for var- 
ious Schwarzchild radii, rg. (b) Crosses give the UV spectral in- 
dices. Curves give the indices predicted by the disk-corona model. 
We see that the (solid) curve for rs = 1.3 X 10 12 cm provides the 
best description of both the UV fluxes and spectra. 



linear than that obtained by extrapolating the X-ray spectra 
to 4 MeV. 

The range of rs obtained above corresponds to M Ss 
8 x 10 6 Af Q . At the largest M and the source size <; 10rs, 
At <; 10 3 s, compatible with the observed longer variabil- 
ity time scales. (Note that M larger than that given above 
can be obtained for a Kerr metric.) The range of the max- 
imum disk temperatures is 7-10 eV. Note that the implied 
range of EF E at 100 keV is 0.1-0.3 keV cm" 2 s" 1 , which is 
much larger than that in 1991-93, when OSSE observed a 
relatively constant 100 keV flux (Z96). 

Summarizing this section, we find that reprocessing by 
a disk is a viable model for NGC 4151 if the X-7 spectra 
are cut off at ~ 100 keV. Our conclusion differs from that 
of PP94, who required that the power-law spectra extend 
to ~ 4 MeV in order to satisfy the energy balance. The 
difference is entirely explained by the different geometrical 
models adopted by PP94 and us. First, our geometry gives 
2(1 + fj,) times more reprocessed flux (for unit observed X-7 
flux) . Physically, this is due to a small solid angle subtended 
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by the disk as seen from a point on the surface of the sphere, 
as well as due to the optically-thick sphere emitting only 
outside (as opposed to the corona emitting both outside and 
inside). Second, our model gives the UV spectrum different 
from that of PP94, with our UV spectrum compatible with 
the data without the need for 7-ray emission harder than 
observed. 

The main caveat for this model is the predicted Comp- 
ton reflection, with R/j, 0.15. Such reflection is only 
marginally allowed by the X-ray data. 



different predictions regarding the specific UV-X-7 variabil- 
ity, and thus can be tested by coordinated multiwavelength 
observations. 

Note that the transmission model produces the UV 
within the X-ray absorber so that the expected absorption 
in the UV is different than that produced in the reflection 
model, where the UV and X-rays both are affected by the 
same absorber (as proposed by Warwick, Smith & Done 
1995). The recent data from the Hopkins Ultraviolet Tele- 
scope appear to rule out the latter possibility (Kriss et al. 
1995), thus favouring the transmission model. 



5 OTHER OBSERVATIONS 

The two models proposed above explain the UV/X-ray cor- 
relation observed in 1983-86. However, the behavior of the 
UV and X-7 in NGC 4151 changes from one epoch to an- 
other. In 1979 May 19-31, EF E (8.5 eV) was ~ 0.25 keV 
cm" 2 s" 1 while EF E (5 keV) was ~ 0.07-0.12 keV cm" 2 s" 1 
(Perola et al. 1982). The UV flux is larger than that implied 
by the specific correlation of Fig. []]. The sampling of the UV 
flux was too sparse to either establish or rule out a UV/X- 
ray correlation in that period with an Fo larger than that in 
1983-86. 

A correlation with a large Fo has been in fact observed 
in 1993 November 30-December 10 (Warwick et al. 1996). 
The EF E (8.5eV) varied within ~ 0.5-0.6 keV cm" 2 s" 1 
[using E(B-V) = 0.05, Section f| together with the 1-2 keV 
flux (measured by ROSAT and ASCA). The EF E (5keV) 
varied in the ~ 0.05-0.25 keV cm" 2 s" 1 range, whereas the 
flux above 50 keV measured simultaneously by OSSE varied 
only within ±15 per cent. 

Preliminary estimates indicate that 1993 December 
data can be explained by either the transmission or reflec- 
tion model; discussion is given in Warwick et al. (1996). We 
only note that a large residual UV flux is required in the 
transmission model whereas the reflection model may ex- 
plain most of the constant UV component as due to disk 
reprocessing of the approximately constant 7-rays. 



6 CONCLUSIONS 

We have found that the X-ray/UV correlation in NGC 4151 
seen by EXOSAT and IV E can be explained by either of 
two models. In the first one, the UV flux is due to reemission 
of the X-ray flux absorbed by small, Thomson-thin clouds, 
which absorption is directly measured in X-rays. The model 
provides an excellent fit to the data. 

In the other model, the X-7 source forms a patchy, dis- 
sipative corona above the surface of an accretion disk. The 
disk reprocesses a fraction of the X-7 flux into the UV, which 
has been proposed before but for a central X-7 source. The 
geometry proposed here allows the UV/X-7 energy balance 
to be satisfied for X-7 spectra cut off at ~ 100 keV, as ob- 
served by OSSE. 

The reflection model requires spectral hardening above 
10 keV, owing to the presence of Compton reflection. This 
can be tested by X-7 observations including coverage of the 
10-50 keV range (e.g., by XTE), where the relative contribu- 
tion from Compton reflection peaks. The two models make 
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